u'"'7  AND  USE  ATJE6MC. 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


UNITED  STATES  AIR  FORCE 
AIR  WEATHER  SERVICE  (NAC) 

AIR  FORCE  GLOBAL  WEATHER  CENTRAL 
OFFUTT  AFB  NE  68113 


y  1 3  ^  7  7 

8°  8  18  J)6 


REVIEW  AND  APPROVAL  STATEMENT 


This  publication  approved  for  public  release.  There  is  no  objection  to 
unlimited  distribution  of  this  document  to  the  public  at  large,  or  by 
the  Defense  Technical  Information  Center  (ETIC)  to  the  National  Technical 
Information  Service (NTIS). 

This  technical  publication  has  been  reviewed  and  is  approved  for  publi¬ 
cation. 

(fL. - 

RONALD  L.  LININGER,  Lt  Colonel,  USAF 
Reviewing  Official 

FOR  THE  CCMtANDER 


ROBERT  C.  SABIN,  Lt  Colonel,  USAF 
Chief,  Technical  Services  Division 


SECURITY  CLASS'FICATION  of  THIS  °fGE  Dote  Entered) 


REPORT  DOCUMENTATION  PAGE 


.  REPORT  SUMOER 

-  80/002 


4.  TITLE  (end  Subtitle) 


flaPM 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  PECl«TKT'-,  CATALOG  NUMBER 


,fA 


5.  TYFE  OF  FEPOF.T  «  PERIOD  COVEREO 


Geomagnetic  Index  Calculation  and  Use  at  AFGWC 


6  performing  org.  report  number 


7.  AUTHORfa) 


8.  CONTRACT  GU  GRANT  NUMBER(t) 


Robert  D.  Prochaska 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

HQ  Air  Force  Global  Weather  Central  (MAC) 
Offutt  AFB,  Nebraska  68113 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

HQ  Air  Force  Global  Weather  Central  (MAC) 
Offutt  AFB,  Nebraska  68113 


^POGRAM  ELEMENT,  PROJECT,  TASK 
I  A  A  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 


April  1980 


13.  NUMBER  OF  PAGES 


4.  MONITORING  AGENCY  NAME  ft  ADDRESS^//  different  from  Controlling  Office)  I  15.  SECURITY  CLASS,  (of  thie  report) 


(same  as  block  9) 


Unclassified 


16.  DISTRIBUTION  STATEMENT  (ol  thl.  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetrmct  entered  in  Block  20.  if  different  from  Report) 


19.  KEY  WOROS  (Continue  on  reveree  eide  If  neceeemry  end  identify  by  block  number) 


Geomagnetic  Index 
Geomagnetic  Field 
Geomagnetic  Storm 
Magnetometer  Analysis 


0.  ABSTRACT  (Conttmre  on  reveree  elde  If  neceeeery  end  identity  by  block  number) 


At  Air  Force  Global  Weather  Central  (AFGWC)  w  use- several  geomagnetic  indices 
to  specify  the  state  of  the  magnetosphe'ra.  T^this  paper  If  introduce*  the  indices 
of  geomagnetic  activity  most  commonly  in  use  by  D0D  agencies,  If discuss  the 
procedures  used  at  AFGWC  to  confute  the  indices  and  briefly  describe  the  uses 
of  each  one. 


DO  |  jan*TS  1473  EDITION  OF  1  NOV  ••  IS  OBSOLETE 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  RAGE  I 


Preface 


x  :r  xr  ~s rs  ‘t**  - 

““*•  -  «•»•  •>  «-•••  *_ 

Branch  compute  several  indices  of  this  activity.  Ifew  members  atZ 

branch  and  staff  meteorologists  at  user  locatiL  noJT 

iar  «ith  geomagnetism.  IM,  rcport  is  .  .  ™  n°m°ny  “* 

for  them.  ^  deslgned  “  “  ^tructionsl  aid 

cal  IrZ^^s  £T  ^T'!ntal  ftctai- 

rt  ~  i 

tcce^er,  W»  Lee  tZ  T  **  *■  «~t  as  of 

-U  not  .cessarilTL  l ^if  o““  **  ClaS’TO™  « 

report  may  be  retained  for  future  10031  Pr°CedUres  chan€e-  This 

lieu  of  operating  instructions.  606006  ^  °0t  **  USed  ** 

Systems  section  of  miCcZlZZZTZ^tT^mZ^  ^ 

r:;  -  - — ^nzrr° 

view  of  the  paper  which  led  to  several  * 

improvements  in  its  content.  several  significant 


Accession  for  a 

HTI6  GRAfcl  W 

DOC  TAB  j 

Unannounced  fl 

Justification 

By_ — — - 

Distribution/., - 

Avail  nbi  \i  tJt~' Cogs3_ 
1  1  Avail  and/or 


TABLE  OP  CONTENTS 


Page 


1.  Introduction  . 1 

2.  Geomagnetic  Field  Strength  Measurements . . . .  2 

a.  Three  Hour,  Single  Station  Indices  . . .  5 

b.  Daily/Flanetary  Indices . . . . .  12 

c.  Other  Geomagnetic  Indices  . . . .. . .  IS 

3*  Use  of  Geomagnetic  Indices  . . 21 

4.  Auroral  Oval  Index,  Qg . 22 

5.  Sumnary . 24 

6.  References  . 26 

7.  Appendix  . . 27 


ii 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  Geodetic  elements  of  the  geomagnetic  field  .  3 

2  Simulated  magnetometer  ouput  .  4 

3  Geomagnetic  Observing  Codes  . . * .  9 

4  UMAGE  code  breakdown  . . 17 

5  Musical  note  diagram  of  Kp . . . . .  19 

6  Gottingen  Ap/Cp  Summary  . . 20 

7  Flowchart  of  AFGWC  geomagnetic  index  calculation  . 25 

LIST  OF  TABLES 

Table  Page 

1  Magnetic  observatories  used  by  AFGWC .  6 

2  Conversion  factors  from  chart  division  to  deflection  ...  6 

3  a  to  K  conversion .  7 

4  Equivalent  amplitude  for  a  given  K  value . .  11 

5  Equivalent  amplitude  conversion  factors  . . 11 

6  Magnetic  observatories  (Gottingen  Ap)  .  13 

7  Minimum  ap  threshold  for  a  given  Kp .  15 

8  Plain  language  Ap  description . 15 

9  Upper  range  limit  in  gammas  defining  Q  values  .  23 

10  Geomagnetic  indices  used  at  AFGWC  23 

iii 


INTRODUCTION 


The  operation  of  several  units  in  the  Department  of  Defense  are  af¬ 
fected  by  geomagnetic  disturbances.  The  effects  of  the  disturbances  in¬ 
clude  increased  atmospheric  drag  on  spacecraft  due  to  heating  of  the 
earth's  atmosphere  and  degraded  high  latitude  communications  due  to  ion¬ 
ospheric  disturbances.  Geomagnetic  observations  are  taken  at  several 
locations.  Air  Force  Global  Weather  Central  personnel  use  them  to  com¬ 
pute  various  indices  of  a  geomagnetic  disturbance  to  meet  the  require¬ 
ments  of  Department  of  Defense  customers. 


1 


GEOMAGNETIC  FIELD  STRENGTH  MEASUREMENTS 


Geomagnetic  field  strength  measurements  are  made  using  an  instrument 
called  a  magnetometer.  The  magnetometer  makes  continuous  measurements 
of  the  instantaneous  geomagnetic  field  strength  along  three  mutually 
perpendicular  axes.  Figure  1  displays  the  conventional  directions  of 
the  axes,  and  the  H,  D,  and  Z  set  or  the  X,  Y,  and  Z  set  of  axes  are 
the  most  common  axes  chosen*.  The  output  of  each  sensor  is  recorded  on 
strip-charts  (see  figure  2).  Geomagnetic  field  strength  is  specified  in 
units  called  gammas.  The  readout  on  the  strip-charts  for  a  sensor  is  in 
units  of  chart  divisions,  where  each  chart  division  represents  a  certain 
number  of  gammas  on  the  axis  of  that  sensor.  At  observing  sites  manned 
by  Air  Weather  Service  (AWS)  personnel,  the  observations  are  made  and 
encoded  in  units  tenths  of  chart  divisions  and  all  work  at  these  observ¬ 
ing  sites  are  in  units  of  chart  divisions.  The  observers  leave  the  con¬ 
version  to  units  of  gammas  to  Air  Force  Global  Weather  Central  (AFGWC). 

The  basic  observing  period  for  geomagnetic  disturbance  measurements 
is  three  hours  long.  The  first  observing  period  of  the  day  begins  at 
00:00:00  universal  time  (UT  or  zulu,  Z)  and  ends  at  02:59:59  UT.  Thus 
there  are  eight  observing  periods  each  day.  A  disturbance  which  occur¬ 
red  at  0143Z  is  included  in  the  0300Z  (00:00:00  UT  to  02:59:59  UT)  ob¬ 
servation.  Observatories  in  the  real-time  network  of  magnetometers  mon¬ 
itored  by  AFGWC  (Table  l)  divide  each  observing  period  into  two  90  min¬ 
ute  sections. 


♦Air  Weather  Service  operated  magnetometers  only  measure  geomagnetic 
field  strengths  along  two  axes,  the  n  and  Z.  It  is  common  practice  in  Air 
Weather  Service  to  refer  to  the  three  sensors  as  "X”,  "Y",  and  Z*  The  "X" 
and  "Yw  sensors  are,  in  reality,  two  sensitivities  of  the  horizontal  component 
measurement.  The  sensors  are  parallel,  with  the  one  called  "Y"  recording  the 
values  at  higher  sensitivity  than  the  "X"  sensor.  The  assumption  made  in  not 
measuring  three  independent  components  is  that  a  disturbance  will  not  affect 
the  geomagnetic  field  strength  only  along  one  axis.  This  assumption,  while 
not  precisely  correct,  is  true  virtually  100$  of  the  time. 
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EARTH'S  CENTER  (+) 


D=DECU NATION  2=  VERTICAL  COMPONENT 

H=  HORIZONTAL  INTENSITY  X=  NORTH  COMPONENT 

Bs TOTAL  INTENSITY  Y=EAST  COMPONENT 

^INCLINATION 


Figure  I.  Geodetic  element*  of  the  geomagnetic  field. 
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Each  magnetometer  observation  gives  the  left  and  right  extent  of  the 
strip-chart  trace  made  during  the  observing  period  (see  figure  2).  Tech¬ 
nically,  this  observation  represents  the  change  of  geomagnetic  field  in^- 
tensity  from  the  quiet  day  curve  (QDC)  along  each  axis,  where  the  QDC  is 
caused  by  currents  due  to  the  apparent  motion  of  the  quiet  sun  (Sq)  and 
moon  (L)  *  In  practice,  of  the  observing  locations  used  at  AFGWC  only 
Boulder  has  a  significant  QDC.  The  analyst  (or  computer  for  automated 
reports)  at  Boulder  removes  the  QDC  prior  to  transmitting  the  observa¬ 
tion;  for  all  other  observing  locations  monitored  by  AFGWC  the  QDC  is  so 
small  that  it  is  ignored. 

AWS  (and  cooperating)  observatories  transmit  their  90  minute  obser¬ 
vations  in  GEMAG  code  (figure  3)*  Each  observation  gives  the  furthest 
left  (maximum  intensity)  and  furthest  right  (minimum  intensity)  value  of 
the  trace  for  each  of  the  three  sensors  during  the  90  minutes.  These 
observations  are  the  basis  for  calculated  indices  of  geomagnetic  distur¬ 
bances. 

Three  Hour,  Single  Station  Indices 

The  basic  index  of  the  disturbance  of  the  geomagnetic  field  is  the 
AMPLITUDE,  a.  This  index  is  the  magnitude  of  the  largest  disturbance 
(from  the  QDC)  of  any  sensor  measured  by  a  magnetometer  during  a  three 
hour  observing  period.  It  is  specified  in  units  of  ganmas  (units  of 
magnetic  intensity) . 

AFGWC  personnel  calculate  an  amplitude  for  each  reporting  period  for 
each  station  listed  in  Table  1.  They  combine  two  consecutive  GEMAG  re¬ 
ports  by  taking  the  larger  of  the  two  maxima  and  the  smaller  of  the  two 
minima.  This  combination  is  the  values  which  would  be  reported  in  a 
three  hour  GEMAG  which  covered  the  same  two  consecutive  90  minute  periods. 
The  difference  between  the  maximum  and  minimum  intensity  for  a  sensor  is 
the  disturbance  level  in  units  of  tenths  of  chart  divisions.  AFGWC  per¬ 
sonnel  convert  this  value  to' the  RANGE,  R,  in  units  of  gammas  (using 
Table  2)  .  The  largest  range  for  the  three  hour  period  is  the  amplitude 
for  that  period. 


Table  1.  Magnetic  observatories  used  by  AFGWC 


Observatory  Station  timber  Geographic  Geomagnetic 


Lat 

Long 

Lat 

Long 

Boulder,  Colorado 

72469 

40  08N 

105  14W 

+49.0 

316.52 

College  Observatory, 
Fairbanks,  Alaska 

25602 

64  52N 

147  50W 

+64.6 

256. 5E 

Goose  Bay,  Labador, 
Canada 

71 816 

55  20N 

60  30W 

+60.5 

11. 9E 

Loring  AFB,  Maine 

72712 

46  57N 

67  53W 

+58.5 

1.5E 

Thule,  Greenland 

04205 

77  29N 

69  10W 

+88.0 

5.7E 

RAF  Upper  Heyford, 

03655 

51  56N 

1  15W 

+50.7 

79.  IE 

England 


Table  2.  Conversion  Factors  from  Chart  Divisions  to  Deflection 


(difference  is 

chart  divisions)  * 

(conversion  factor) 

=  (range  in  gammas) 

Observatory 

1st  Sensor 

2nd  Sensor 

3rd  Sensor 

Boulder 

10.0 

10.0 

* 

College 

80. 5 

80. 4 

81.5 

Goose  Bay 

125.0 

20.0 

100.0 

Loring 

125.0 

12.5 

50.0 

Thule 

50.0 

10.0 

20.0 

Upper  Heyford 

50.0 

12.5 

50.0 

*  sensor  not  used/not  available 
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GEMAG 

IIIII  IJWDD  3GGgg  4S/NN 

11111  20MN1  30MX1  40MN2  50MX2  60MN3  70MX3  99999 

22222  20MN1  30MX1  40MN2  50MX2  60MN3  70MX3  99999 

IIIII  ■  station  number 

UfflDD  =  year/month/day;  e.g.,  1  Apr  1979  is  90401 

GGgg  =  hour/minute  of  the  end  of  the  observation  period 

s  =  report  status  (9  =  data  valid  for  90  minutes  ending  at  GGgg) 

11111  =  this  line  is  for  90  minutes  ending  at  GGgg 

22222  =  this  line  is  for  90  minutes  just  prior  to  11111  line  and  is 

optional.  If  used,  omit  the  99999  group  at  the  end  of  the 
previous  ( 11111 )  line. 

MN1  =  Minimum  pen  deflection  on  sensor  1 

MX1  =  Maximum  pen  deflection  on  sensor  1 

MN2/MX2  m  minimum/amximum  for  sensor  2 

MN3/MX3  =  minimum/maximum  for  sensor  3 

a.  GEMAG  Code 


MAGTR  IIIII  DDHH/  GGGll  99999 

I Till  »  station  number 

DDHH  m  day/hour  of  start  of  3  hour  period 

GGG  =  amplitude,  a 

NOTE:  K  index  is  normally  appended  to  this  message  as  a  plain  language 
comment • 

b.  MAGTR  Code 


Figure  3.  Geomagnetic  Observing  Codes.  These  codes  are  used  for  real- 
time  reports  of  geomagnetic  observations. 
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Amplitude  is  a  linear  index  of  geomagnetic  disturbance  level.  An 
amplitude  of  50  represents  a  geomagnetic  disturbance  (in  gammas)  twice 
as  severe  as  that  represented  by  a  25* 

Amplitude  is  strongly  dependent  on  the  observatory  location  in  geo¬ 
magnetic  coordinates.  Instruments  located  in  auroral  latitudes  react 
most  strongly,  and  those  near  30°  north  or  south  react  least  strongly  to 
a  disturbance  of  the  worldwide  geomagnetic  field.  You  must  remove  the 
geomagnetic  latitude  dependence  of  the  amplitudes  before  you  can  compare 
observations  from  different  observatories. 

In  1938,  the  K  scale  was  adopted  for  the  observatory  at  Niemegk 
(near  Potsdam,  GDR)  as  a  quasi-logarithmic  index  of  geomagnetic  activity. 
The  values  chosen  gave  a  reasonable  distribution  of  frequency  of  occur¬ 
rence  across  the  scale  from  zero  (very  quiet)  to  nine  (violently  distur¬ 
bed)  •  All  observatories  could  not  use  the  same  scale  to  convert  ampli¬ 
tude  to  Kf  the  "LOGARITHMIC”  AMPLITUDE;  otherwise  auroral  stations  would 
show  much  greater  K  values  than  would  lower  latitude  stations  for  the 
same  storm.  So  to  get  a  conversion  from  a  to  K  for  some  stations  we 
remove  the  geomagnetic  latitude  dependence.  The  1938  data  for  any  pai>- 
ticular  station  was  compared  to  that  from  Niemegk  to  get  the  same  fre¬ 
quency  of  occurrance  of  K  in  each  category  at  each  station.  As  new  ob¬ 
servatories  are  established  their  amplitude  values  are  compared  to  cal¬ 
ibrated  observatories  [7]*  Table  3  is  used  to  convert  a  to  K  for  the 
observatories  which  report  to  AFGWC  in  real-time. 

Since  the  K-index  is  based  on  a  single  station,  it  represents  region¬ 
al  (longitude)  conditions  and  will  include  local  variations  in  geomag¬ 
netic  activity,  including  the  Sq  and  L  variations.  A  possible  index  to 
express  the  geomagnetic  disturbance  on  a  worldwide  scale  might  be  a  sim¬ 
ple  average  of  the  K  indices  from  several  (or  all  available)  stations. 
However,  the  K  index  is  not  linear,  so  a  simple  average  would  not  be 
useful  even  if  we  had  a  uniform  geographic  distribution  of  observatories. 
The  first  step  to  a  planetary  index  is  to  compute  a  single  station 
"equivalent"  index. 
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The  EQUIVALENT  AMPLITUIB,  a^,  is  a  1  .near,  single— station,  three  hour 
index  which  has  had  the  geomagnetic  latitude  dependence  removed.  The  in¬ 
ternational  way  of  determining  a^  is  to  use  Table  4,  converting  the  ampli¬ 
tude  to  K  and  K  to  the  equivalent  amplitude.  AT  AFGWC,  the  analyst  con¬ 
verts  amplitude  directly  to  equivalent  amplitude,  dividing  each  stations 
amplitude  by  a  latitude  dependent  factor  (see  Table  5)*  Since  a^  is  lin¬ 
ear  and  station-independent  9  various  values  can  be  averaged  for  several 
time  periods  to  get  other  indices;  e.g.f  the  a^  values  during  a  geomagnetic 
disturbance  can  be  averaged  to  determine  the  severity  of  the  disturbance. 

EXERCISE  1.  Using  Tables  1-3  and  5y  figure  3y  said  the  following  reports: 


GEMAG 

71316 

11111 

90403 

20449 

31030 

30403 

49/12 

40497 

50215 

60406 

70337 

22222 

20452 

30444 

40528 

50474 

60335 

70328 

99999 

GEMAG 

71816 

11111 

90403 

20/// 

31200 

30/// 

49/06 

404&& 

50467 

60330 

70329 

99999 

complete  the  following: 
station  name: 

date  of  reports:  . 

0730Z  -  1030Z  amplitude  =  _ 

0730Z  -  1030Z  K  index  =  _ 

0900Z  -  1200Z  amplitude  =  _ _______ 

0900Z  -  1200Z  K  index  =  _ 

0900Z  -  1200Z  equivalent  amplitude  =  _ 

EXERCISE  2.  MAGTR  is  a  code  which  is  used  by  Boulder  to  report  a/K  val¬ 
ues  directly.  Use  figure  3  to  interpret  the  following  report. 

MAGTR  72469  1215/  22711  99999 
station:  Boulder  Day  ___ ___ 

Time:  Z  to  Z,  a  «  t  K  »  _ 
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Table  4*  Equivalent  Amplitude  for  a  Given  K  Value 


KOI 
aR  0  3 


2  3  4  5 

7  15  27  48 


6  7  8  9 

80  140  240  400 


Table  5*  Equivalent  Amplitude  Conversion  Factors 
(a)  -  (conversion  factor)  =  aj( 


Observatory 

Boulder 
College 
Goose  Bay 
Loring 
Thule 

Upper  Heyford 


Factor  (gammas) 

2.0 

10.0 

6.5 

3.0 

4.0 

2.4 
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Daily/Planetary  Indices 

The  indices  which  we  have  examined  this  far  are  valid  only  for  one 
station  and  they  are  valid  during  a  particular  three  hour  period.  We 
will  now  examine  indices  used  to  specify  geomagnetic  activity  level  on  a 
daily  and/or  planetary  scale. 

The  EQUIVALENT  DAILY  AMPLITUIE,  A&,  specifies  the  geomagnetic  activ¬ 
ity  level  for  a  particular  station  on  a  particular  day.  A^  is  obtained 
by  averaging  the  chosen  station’s  eight  values  for  that  day.  Since 
aft  is  independent  of  the  station’s  geomagnetic  latitude,  can  be  used 
to  compare  the  daily  geomagnetic  activity  at  various  observatories. 

When  discussing  A^  values  for  various  observatories,  it  is  common  prac¬ 
tice  to  replace  the  subscript  K  with  a  two  letter  subcript  to  identify 
the  station.  For  example,  the  %  for  Fredricksburg,  Virginia  is  noiv- 
mally  denoted  as  App  (called  A-FEED  by  NOAA  personnel  at  Boulder)  • 

Since  AK  is  linear,  we  can  average  several  days  worth  to  get  weekly, 
monthly,  or  even  yearly  values  of  the  average  disturbance  at  that  sta¬ 
tion.  Any  one  station  index  does  not  specify  the  worldwide  level  of 
distribution  of  the  geomagnetic  field,  so  we  will  next  examine  the  so- 
called  planetary  indices,  indices  of  worldwide  geomagnetic  activity. 

All  planetary  indices  of  geomagnetic  activity  are  computed  by  aver¬ 
aging  selected  single  station  indices.  Ideally,  one  should  use  a  large 
number  of  stations  which  are  evenly  distributed  about  the  globe.  In 
practice,  a  small  number  of  observations  taken  at  selected  locations  are 
used.  By  averaging  the  aK  values  from  twelve  specific  stations  (see 
Table  6) ,  the  personnel  of  the  International  Service  for  Geomagnetic 
Indices  determine  the  EQUIVALENT  PLANET  ART  AMPLITUDE,  ap.  This  value 
represents  the  official  value  of  the  activity  of  the  entire  geomagnetic 
field  for  the  three  hours.  However,  this  index  is  not  available  in  real 
time  and  so  is  not  routinely  used  by  AFGWC#  Personnel  at  AFGWC  normally 
refer  to  it  as  the  Gottingen  ap. 


12 


gj 

On 

bO 

C  <0 
Q  LA 
O  >3  CA 

£ 

§ 

fe  8 

£  la  o 

O  M  VO 


s 

w 

w 

££) 

w 

« 

Cx} 

W 

w 

W 

o 

LA 

3 

Sv 

3 

o 

3 

rH 

CA 

3 

to 

CA 

3 

C- 

rH 

ct 

CV 

CV 

rH 

« 

CA 

a 

On 

NO 

LA 

LA 

CA 

§ 

cv 

8 

CV 

rH 

IS 

rH 

3 

55 

rH 

8 

525 

F- 

»5 

LA 

R 

g 

R 

R 

8 

cv 

O 

IA 

rH 

CA 

-4“ 

-4 

o 

-4 

rH 

o 

ON 

fN- 

LA 

LA 

-4 

CA 

cv 

rH 

CA 

to 

CA 

LA 

LA 

VA 

LA 

LA 

LA 

LA 

LA 

-4 

CA 

IS 

e 

0) 


■3 


rH 

m 


rH 

-P 

I 

CO 


■s 


I  3 


eO 

«  i 


i 


a  a 


I 

§ 

a 


o 

cO 


cd 

■§ 


1  I 


■g  <? 

- 1 


! 

8 

IS 

£ 


U 

JS 

-P 

£ 


of 

O 

8 

3 

0> 

i 

•* 

H 

<D 

«P 

<§ 

S 

■P 

H-> 

•rl 

CO 

a 

£ 

£ 

I 

$ 

£ 

<§ 

£ 

£ 

£ 

$ 

to  no 

<J  8 

•H 

-P 

fl) 


to 


LA 

o 


-4 

LT\ 

R? 


E  vr\  rH  o 

O  -P  •  •  • 

^  tfl  CM  tO  O 

u  *3  ja  ^ 


LA  O 


to 

ON 


Si 

to 


o 

(A 


2 


ON 


O  O  On  LA 

•  *  •  • 

ON  t>-  q>  CV 

C^-  -4  -4  la 

CA  CA  C\J 


to 

LA 

to 

LA 

CV 

NO 

o 

• 

• 

• 

• 

• 

• 

• 

LA 

to 

rH 

-4 

-4 

-4 

LA 

LA 

LA 

sp 

LA 

LA 

LA 

LA 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

NO 


f 


13 


AFGWC  personnel  use  a  network  of  real  time  stations  (Table  l)  to  sim¬ 
ulate  the  EQUIVALENT  PLANETARY  AMPLITUEE,  ap •  They  average  the  a^  values 
for  five  of  these  stations  (excluding  Thule)  to  get  a  value  which  very 
closely  approximates  ap  and  which  at  AFGWC  is  called  ap.  In  the  remain¬ 
der  of  this  discussion  I  will  refer  to  the  AFGWC  value  as  ap,  although 
my  comments  are  more  precisely  correct  for  the  Gdttingen  value-  ap,  a 
linear  index,  ranges  from  zero  to  well  over  one  hundred. 

Thule,  like  other  polar  cap  stations,  is  excluded  from  planetary 
calculations  because  it  has  a  poor  correlation  to  worldwide  disturbances. 
Polar  stations  respond  very  strongly  to  disturbances  of  the  geomagnetic 
field  which  affect  only  the  polar  regions.  They  respond  less  strongly 
to  worldwide  disturbances  than  do  auroral  stations.  Thus,  they  are  ex¬ 
cluded  from  planetary  calculations,  as  they  would  bias  the  calculated 
index  away  from  the  "planetary"  value. 

The  "LOGARITHMIC"  EQUIVALENT  PLANETARY  AMPLITUEE,  Kp,  is  a  three 
hour  value  of  geomagnetic  activity  for  the  entire  earth,  using  a  quasi- 
logarithmic  scale.  Kp  values  range  from  0  to  9  in  28  steps;  i.e.,  0,  0+, 
1-,  1,  1+, ...7+,  8-,  8,  8+,  9-*  9*  At  AFGWC,  the  conversion  from  ap  to 
Kp  is  made  using  Table  7* 

The  specification  of  the  total  planetary  level  of  geomagnetic  acti¬ 
vity  during  a  24-hour  period  is  an  average  of  the  eight  ap  values  for 
that  period.  The  resultant  index  is  the  DAILY  EQUIVALENT  PLANETARY  AM¬ 
PLITUDE,  Ap.  At  AFGWC  the  classification  of  geomagnetic  disturbances  is 
based  on  the  Ap  and  the  plain  language  descriptions  listed  in  Table  8. 
Again,  as  with  ap,  the  official  value  is  the  GSttingen  Ap,  calculated 
from  the  stations  listed  in  Table  6,  while  AFGWC  calculates  a  similar 
value  using  the  real  time  reporting  stations. 
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Table  7.  Minimum  ap  Threshold  for  a  Given  Kp  value 


aP 

Kp 

*P 

Kp 

ap 

Kp 

0 

0 

15 

3 

80 

6 

2 

0+ 

18 

3+ 

94 

6+ 

3 

1- 

22 

V* 

111 

7- 

4 

1 

27 

4 

132 

7 

5 

1+ 

32 

4+ 

154 

7+ 

6 

2- 

39 

5- 

179 

8- 

7 

2 

5 

207 

8 

9 

2+ 

56 

5+ 

236 

8+ 

12 

3- 

67 

6- 

300 

9- 

400 

9 

Table  8.  Plain  Language  Ap  Description 

Condition  Minimum  Ap 


Quiet  0 
Unsettled  7 
Active  15 
Minor  Storm  30 
Major  Storm  50 
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EXERCISE  3*  Given  the  following  amplitudes  (in  ganmas) 


End  Time 

Station 

02/00 

02/03 

02/06 

02/09 

02/12 

02/15 

02/18 

02/21 

02/24 

Loring 

1 

0 

12 

3 

90 

870 

210 

90 

180 

College 

15 

20 

50 

30 

550 

4000 

2300 

650 

1100 

Goose  Bay 

13 

13 

33 

52 

358 

2275 

1073 

455 

942 

Boulder 

0 

0 

10 

4 

70 

520 

90 

50 

100 

Upper  Heyford  3  1 

Complete  the  following: 

14 

3 

60 

528 

156 

120 

204 

End  Time 

02/00 

02/03 

02/06 

02/09 

02/12 

02/15 

02/18 

02/21 

02/24 

Station 

Loring 

College  aj^ 

Goose  Bay  aj( 

Boulder  a^ 

Upper 

Heyford 

ap 

Kp 

What  is  the  %  for  College  for  day  02?  For  Goose  Bay? 

What  is  the  Ap  for  the  day  (02)? 

EXERCISE  4*  Encode  the  Boulder  MAGTR  for  the  time  period  ending  02/1900Z, 
using  the  data  from  Exercise  3 • 

EXERCISE  5*  URSI  (international  Radio  Science  Union)  observatories  trans¬ 
mit  daily  summaries  of  their  K  index  values  and  the  Ajj  value  for  that  day 
using  the  code  UMAGE  (figure  4)*  Use  the  data  from  exercise  3  to  encode  a 
UMAGE  for  Loring  from  01/21Z  to  02/21Z,  and  one  for  Upper  Heyford  from 
02/00Z  to  02/24Z  (indicates  SSC  of  02/0315Z  for  both  stations)*  Since 
Loring  and  Upper  Heyford  are  AWS  stations*  they  do  not  calculate  the  in¬ 
dices  included  in  the  UMAGE  reports* 
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UMAGE  IIIII  JJTTa  bbkkk  kkkkk  cHHran  dHHmm  eeeee 


IIIII  =  station  number 

JJTT  -  day/hour  of  start  of  24  hour  period  reported 

a  a  check  sum  last  digit  (this  is  the  last  digit  of  the  sum  of  the  2-b 
and  8~k  values) 

bb  =  last  two  digits  of  Ag  for  day  JJ  (use  //  or  XX  if  not  available) 
k  =  K  index,  one  for  each  3-hour  period  beginning  at  JJTT 
c  =  special  phenomenon 

1  =  storm  end 

2  *  bay  (PSC) 

3  =  typical  crochet 

4  =  provisional  figures 

6  =  gradual  storm  beginning 

7  =  SSC 

8  =  very  marked  SSC 

HHmm  a  time  of  special  phenomenon 
Optional  Groups 

d  -  5  means  special  phenomenon  of  H-  component  minimum 
eeeee  =  intensity  of  H-component  minimum  in  gammas 


Figure  4*  UMAGE  Code  Breakdown 
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Other  Geomagnetic  Indices 


There  are  several  other  indices  of  geomagnetic  activity.  Indeed, 
each  geomagnetic  index  (including  those  previously  discussed)  was  de¬ 
veloped  for  a  particular  use,  and  the  relationship  between  indices  is 
generally  not  a  simple  one.  The  people  who  use  a  particular  index  of¬ 
ten  have  a  strong  bias  toward  that  index  as  "the"  indicator  of  geomag¬ 
netic  activity;  I  will  avoid  any  indorsement  of  a  particular  index.  I 
will  now  briefly  introduce  the  more  generally  accepted  geomagnetic  in¬ 
dices  which  are  not  commonly  used  by  DOD  agencies. 

The  Kp  index  may  be  obtained  from  K  without  the  computation  of  aj(. 

By  removing  the  longitudinal,  or  local  time,  variation  effect  and  the 
seasonal  effect  from  a  K  index,  you  can  obtain  the  STANDARDIZED  K  index, 
Kg.  Kp  is  the  average  of  the  Kg  values  from  the  12  stations  used  to 
compute  Ap*  Kp,  by  the  way,  is  due  to  a  combination  of  high  and  low 
latitude  effects  and  the  correlation  between  it  and  some  other  geomag¬ 
netic  indices  (like  AE  and  Dst  below)  is  not  a  simple  relationship.  Kp 
is  published  in  the  so-called  musical-note  format  shown  in  Figure  5» 

The  magnetic  character  Figure,  C,  is  a  subjectively  obtained  index 
for  a  24  hour  interval.  Each  participating  observatory  evaluates  its 
daily  magnetograms  and  assigns  one  of  three  values:  zero  for  quiet, 
one  for  moderately  disturbed,  or  two  for  very  disturbed.  An  average  of 
the  contributing  stations,  to  the  nearest  tenth  of  a  value,  is  called 
Ci*  Variations  of  this  index  include  Cp,  which  is  based  on  conversion 
from  Ap,  and  C9,  which  ranges  from  0  to  9*  Ap  and  Cp  are  published  in 
the  format  of  Figure  6. 

Quantitative  studies  use  indices  which  measure  the  disturbance  by 
region.  Dst  and  AE  are  the  most  conmonly  used  of  these  indices.  The 
storm  time  disturbance  index,  Dst,  is  the  longitudinal  average  depression 
of  the  low  latitude  horizontal  component  in  units  of  ganmas.  This  index 
is  approximately  proportional  to  the  total  kinetic  energy  of  the  parti¬ 
cles  which  are  injected  and  trapped  in  the  radiation  belts  during  the 
main  phase  of  the  geomagnetic  storm*  The  AURORAL  ELECTROJECT  index,  AE, 
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DAYS  IN  SOLAR  ROTATION  INTERVAL 


Figure  5«  Musical  note  diagram  of  Kp 


IUGG:  ASSOCIATION  OF  GEOMAGNETISM  AND  AERONOMY 
(International  Service  of  Geomagnetic  Indices) 

International  quiet  and  disturbed  days? 

Quietest  Days  1  -5:  17  16  6  31  3 

Quietest  Days:6-10:  10A  4A  5A  8A  13A 

Most  Disturbed  Days  1-5:  22  25  24  26  19 

GEOMAGNETIC  PLANETARY  INDICES 
Three-hourly t  Kpj  Daily :  Ap  and  Cp 


May  1979 
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Figure  6.  Gottingen  Ap/Cp  Summary 
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is  a  monitor  of  the  intensity  of  the  auroral  electroject,  based  on  hori¬ 
zontal  magnetic  component  measurements  taken  by  a  set  of  auroral  zone 
magnetic  stations.  The  set  of  stations  should  be  uniformly  disturbed  in 
longitude  to  cancel  the  local  time  dependence  of  the  index-  Dst  is  typi¬ 
cally  specified  With  hourly  time  resolution  while  AE  normally  is  specified 
with  two  to  five  minute  resolution,  [l] 

USE  OF  GEOMAGNETIC  INDICES 

The  three  hour,  single  station  indices  (a,  K,  aj^)  were  designed  for 
use  in  classifying  short-lived  activity  in  a  small  area.  The  amplitude, 
because  of  its  dependence  on  geomagnetic  latitude,  is  of  use  only  in  com¬ 
paring  various  observation  periods  at  the  same  station.  The  equivalent 
amplitude  a^,  and  the  K-index  are  used  to  compare  observations  at  various 
locations,  as  the  geomagnetic  latitude  dependence  of  each  has  been  re¬ 
moved.  The  K-index  ranges  only  from  zero  to  nine  while  the  value  can 
reach  a  few  hundred;  the  K-index  is  more  commonly  used  than  aK-  Small 
(spatial)  scale  disturbances  of  the  geomagnetic  field  which  typically 
last  one  to  four  hours  are  called  substorms;  the  K-index  is  widely  used 
when  substorms  are  examined. 

The  equivalent  daily  amplitude,  %  is  used  to  classify  the  activity 
level  on  various  days  at  a  single  station.  Since  eight  3-hour  a^  values 
are  averaged,  the  effect  of  an  individual  substorm  is  minimized,  and  the 
primary  factor  classified  by  A&  is  the  overall  character  of  the  geomag¬ 
netic  field  at  that  location*  Using  the  appropriate  A%  indices,  we  can 
compare  the  severity  of  a  particular  geomagnetic  storm  as  seen  at  var¬ 
ious  locations,  or  compare  the  geomagnetic  activity  of  various  days  as 
observed  at  a  single  location. 

We  use  multi-station  indices  to  specify  the  overall  state  of  the 
magnetosphere.  The  equivalent  planetary  amplitude,  ap,  and  its  loga¬ 
rithmic  version,  Kp,  define  which  phase  of  a  geomagnetic  storm  (a  major 
disturbance  of  the  magnetosphere)  is  in  progress  during  a  particular 
time  period.  Typically  the  largest  Kp  values,  occur  during  the  main 
phase  of  a  storm,  while  slowly  declining  Kp  values  occur  during  the 
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recovery  phase.  The  daily  equivalent  planetary  amplitude,  Ap,  is  used  to 
classify  the  overall  intensity  of  a  geomagnetic  storm  (see  Table  8). 

AURORAL  OVAL  INEEX  Qe 

The  index  Q  is  a  measure  of  the  aize  of  the  polar  cap  and  the  auroral 
oval  around  it.  Q,  the  polar  range  index,  varies  with  the  intensity  of 
the  ionospheric  currents  which  flow  across  the  polar  cap.  Stronger  polar 
cap  currents  will  produce  a  larger  auroral  oval.  Q  is  a  quasi-logarith- 
mic  index  which  can  be  computed  only  for  observatories  located  poleward 
of  5 8°  geomagnetic  latitude  (see  Table  9)*  Q  is  a  15  minute  index  based 
on  variations  of  the  horizontal  magnetic  field  strength  (vector  H  in 
Figure  l) .  It  is  a  short-term  index  and  reacts  very  strongly  to  geomag¬ 
netic  substorms. 

AFGWC  personnel  do  not  compute  Q,  but  use  an  equivalent  Q  index,  Qg» 
to  monitor  the  auroral  oval  size  and  location.  They  analyze  the  equators 
ward  boundary  of  the  diffuse  aurora  seen  on  DMSP  imagery,  specifying  the 
location  and  time  of  occurrence.  Computerized  analysis  of  precipitating 
electron  measurements  taken  by  EMSP  also  results  in  the  location  and 
time  of  auroral  occurrences  due  to  those  particles.  These  times  and  lo¬ 
cations,  in  Corrected  Geomagnetic  coordinates  [10]  are  used  to  specify 
the  Q  value  which  would  have  caused  that  auroral  oval  extent.  AFGWC 
personnel  use  the  Qg  values  with  Whalen’s  auroral  oval  plotter  [10]  to 
monitor  auroral  and  polar  disturbances. 


22 


Table  9. 

Upper 

Range 

Limit 

in 

Gammas 

Defining  Q  Values 

Q 

0  1 

2 

3 

4 

5 

6  7  8 

9 

10 

Upper 

Limit 

10  20 

40 

80 

140 

240 

400  660  1000 

1500 

2200 

Table 

10.  Geomagnetic 

Indices  used  at 

AFGWC 

Index 

Time  Scale 

Spatial  Scale 

Comnents 

a 

3  hour 

local 

K 

3  hour 

local 

Pseudo-logarithmic 

aK 

3  hour 

local 

Geomagnetic  Latitude  Independent 

ak 

24  hour 

local 

Geomagnetic  Latitude  Independent 

aP 

3  hour 

planetary- 

KP 

3  hour 

planetary 

Pseudo-logarithmic 

AP 

24  hour 

planetary 

Q 

15  minute 

polar  cap 

Not  used  directly,  but  inferred 
from  auroral  observations 
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SUMMARY 


Geomagnetic  activity  indices  are  measures  of  the  level  of  disturbance 
of  the  geomagnetic  field.  Magnetometer  observations  are  local  measure¬ 
ments  of  geomagnetic  disturbances.  AFGVTC  personnel  use  the  measurements 
from  a  network  of  observing  locations  to  calculate  the  indices  Kf  ax»  apf 
Kp  and  Ap  (see  figure  7).  DoD  personnel  use  these  and  the  indices  %  and 
Qg  to  monitor  terrestial  magnetospheric  activity.  Comparison  of  appro¬ 
priate  indices  reveals  the  short  and  long  term,  local  and  worldwide  scale 
and  severity  of  a  geomagnetic  disturbance  (see  Table  10). 
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APPENDIX 


Answers  to  Exercises 
Exercise  1 

Station  71 8l  6  -  Goose  Bay  Labrador 


Date  of  Report  -  90503  -  3  April  1979 


TIME  PERIOD 

X-MAX  X-MIN 

Y-MAX 

Y-MLN 

Z-MAX 

Z-MIN 

0730-0900Z 

44.9  40.3 

49.7 

21.5 

40.6 

33.7 

09001030Z 

45.2  44*4 

52.8 

47.8 

33.5 

32.8 

1030-1200Z 

n/r  n/r 

48.8 

46.7 

33.0 

32.8 

N/R  indicates  none  reported 

TIME  PERIOD 

Range  (Chart  Divisions) 

Range  (Gammas) 

X  Y 

Z 

X 

Y 

Z 

0730-1030Z 

4.9  31.3 

7.8 

612.5 

626.0 

7 80.0 

0900-1200Z 

*  6.1 

0.7 

« 

122.0 

70.0 

TIME  PERIOD 

Amplitude 

K 

Equivalent  Amplitude 

0730-10302 

7  80 

7 

120 

0900-12002 

122 

3 

18.8 

Exercise  2 


Day  12th  of  Month 
Time  period  1500Z  to  1800Z 
a  »  227 
K  -  7 
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Exercise  3 

(values  to  nearest  whole  number) 

End  Time 

02/00 

02/03 

02/06 

02/09 

02/12 

02/15 

02/18 

02/21 

02/24 

Station 

Loring 

0 

0 

4 

1 

30 

290 

70 

30 

60 

College 

1 

2 

5 

3 

55 

400 

230 

65 

110 

Goose  Bay 

2 

2 

5 

8 

55 

350 

165 

70 

145 

Boulder 

0 

0 

5 

2 

35 

260 

45 

25 

50 

Upper 

Heyford 

1 

0 

6 

1 

25 

220 

65 

50 

85 

ap 

1 

1 

5 

3 

40 

304 

115 

48 

90 

KP 

0 

0 

1+ 

1- 

5- 

9- 

7- 

5 

6 

A  College  *  ^9  ^GooBe  Bay 

=  100 

Ap  =  76 

Exercise  4 

MAGTR  72469 

0215/ 

04511 

99999 

Exercise  5 

UMAGE 

72712 

01214 

XX001 

04964 

70315 

UMAGE 

03655 

02002 

56020 

48656 

70315 
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